Analysis of Neutrino Signals from SN1987A 
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Qs . Summary. We study SN1987A neutrino events through a likelihood analysis 

r\J ' with one-component (cooling) and two-component (accretion and cooling) emission 

model. We show that there is a 3.2a hint for the initial accretion phase. 

' ■ 1 Introduction 

o 

H 

"^ ' On 23'"'* February 1987 four neutrino detectors collected a burst of events from 

C3 . supernova (SN) explosion. The signal was detected by Kamiokande-II (KII) 

in Japan (11,-1-5 below threshold, events) [1]; 1MB in Michigan (8 events) [2] 
and Baksan in Russia (5 events) [3], for an amount of 29 neutrino events in a 
^ ' window of T = 30 sec. Four hours before 5 other events in a time window of 

OO , 7 sec were detected by the LSD experiment [4]. To the best of our knowledge, 

the first phase of the neutrino emission, revealed by LSD, cannot be explained 
inside the standard description of a core-collapse SN [5] , whereas the second 
main phase of the neutrino emission can be described using the standard 
scenario that we adopt in the present work. 

QQ ' We investigate the neutrino events detected by KII, 1MB and Baksan, in 

f^ , order to obtain maximal information through a detailed statistical analysis. 

We reconstruct the likelihood function for all events and maximize this proba- 
bility varying some theoretical parameters related to the emission models that 

^^ ' we take into account. In our analysis, we consider the time-energy distribution 

of the signal, the directional information of the data, the detectors properties 
(e.g. the efficiency functions). Furthermore we include a detailed description 
of the background, following and improving the seminal work of Loredo and 
Lamb [6, 7]. 
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2 Emission models 



To grasp the emission models used in our data analysis, now we briefly describe 
the main phases of the so-called standard model for a core-collapse SN [8]. 
At the end of its life, a massive star consists of a sequence of concentric 
shells corresponding to the relics of different burning phases and its inner 
core is formed by iron, which is the final stage of nuclear fusions. The iron 
core grows due to silicon shell burning and, when it reaches a mass of about 
IAAMq (the Chandrasekhar limit mass), the electron degeneracy pressure 
can not support the structure's weight and then the core collapses. At the 
densities and temperatures involved, the processes of electron capture, (3 decay 
and partial photodisintegration of iron-group nuclei to alpha particles occur 
and cause the acceleration of the collapse and the neutrino trapping in the 
core. This collapse proceeds until nuclear densities of about 10^'^g/cm^ are 
reached. At this density the nuclear matter is nearly incompressible and the 
outer iron core rebounds driving a shock wave that propagates outwards, 
whereas the external region falls inwards at supersonic speed. The explosion 
mechanism is still uncertain but, in 'delayed scenario' [9], the shock seems 
to loose its energy because of the dissociation of heavy nuclei into nucleons 
and because of the neutrino emission, that grows when the shock crosses the 
neutrinosphere. The weakened shock stagnates and transforms into a standing 
accretion shock, whereas the outside matter falls inward and joins the nascent 
compact remnant. We call this phase accretion (suffix a) phase and we suppose 
that, in this phase, Ve and Ve are produced in similar amount by the ep —^ nve 
and e^n -^ pv^ processes. The accretion phase occurs within the first second 
of Ve emission, and for this non-thermal phase [10] wc assume the following 
parameterized neutrino flux 
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where Yn — 0.6, Ma is the accreting mass exposed to the positrons thermal 
flux, g{E,T) ~ £^2/ [1 + exp(i?/r)], a Fermi-Dirac distribution with tem- 
perature Tq, ag^+n{Ev) is the cross section of positron interactions increasing 
quadratically with E^. The time scale of accretion process (namely To) ap- 
pears in the following function e{t) = exp[— (i/To)^°]/[l -|- i/(0.5 s)]. 
Taking into account the neutrino oscillations and the accretion assumptions 
(i.e. '^^accr,u = 0)j the D^, accretion flux is reduced by the factor Pee = cos^Oi^ 
(the survival probability function [11, 12]), hence the total flux becomes 



<P (f E ) — P ■ <P^ (f E 



(2) 



The nascent proto- neutron star evolves in a neutron star (with radius Rns) 
and this process is characterized by an intense flux of all the species of 
(anti)neutrinos. We call this phase cooling phase (suffix c), a thermal phase 
with a longer time scale, and we suppose that an equal amount of energy goes 
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in each species (equipartition hypothesis). The adopted form of the parame- 
terized antineutrino flux, differential in the energy is 

1 TTC 

'^"ooiit, E,) - i^^^^ [^^^l 9{E., nt))] , (3) 

that is a standard black body emission from the neutrinosphere with radius 
Re ~ Rns- The time scale of cooling emission Tc is included in the function 
T{t) = Tc exp[— i/(4Tc)]. In this phase all species of neutrino are emitted and 
we have to consider the neutrino oscillations to obtain the total Ve flux at the 
detectors. Assuming equipartition and normal hierarchy for neutrino masses, 
the total cooling flux of the electron antineutrino is 

^coolit, E,) = Pee<o;,P. {t, E,) + (1 - Pee)€ooL,.. (*, E,). (4) 



3 Data analysis 

Let us construct the likelihood function for the data set. We assume that the 
detected v^ interact for Inverse Beta Decay (IBD), i>ep -^ ne^ . The signal 
rate (triply differential in time, in the positron energy Ee, in the cosine of the 
angle between the antineutrino and the positron directions) is given by 



^-^{E.,co^e) 7^d{Ee) Ucose) $^At,E,) ^ 

a cos a abjp 



S{t,Ee,cose)^Np .J— ^{E,, cos 6) 7^d{Ee) UcosO) $^At,E,) -^, (5) 



where Np is the number of targets (free protons) in the detectors, a is the IBD 
cross section [13], r]d is the -detector dependent-average detection efficiency, S,d 
is the angular bias {S,d — 1 for Kamiokande-II and Baksan whereas ^^(cos 9) = 
1+0.1 cosO for 1MB [14]) and, finally, ^p_. is the total fiux of i^e, sum of the two 
term (Pacc and 'Pcooi shown respectively in Eqs. (2) and (4). The theoretical 
parameters, that we have to deduce by fitting the data, are included in this 
last term and are 6 parameters (3 for each phase): Ma, Ta and Tq for accretion; 
Re, Tc and Tc for cooling. The likelihood function is 

E= U ^d, (6) 

d—k.i,b 

where the suffix d ranges over the detectors. Using the Poisson statistic, the 
hkehhood function for each detector is 

^d - . -*/(*+*^''^* J|e^(*.+*.)- ^+ls{U+td,Ec,c,)G.{Ec)dEc 

(7) 
where Nd is the number of events for each detector and the suffix i refers to 
the i-th event {i = l...Nd)- The time td, called "offset time", is the temporal 
gap between the arrival of the first neutrino to the Earth and the detection of 
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first neutrino event in the detector. As consequence, we add 3 new parameters 
td to find out by our data analysis. The term /jj is the detector live fraction 
{fd = 1 for KII and Baksan, whereas fd = 0.9055 for 1MB), r^ is the detector 
dead time (r^ = for KII and 1MB, whereas Td = 0.035 for 1MB). Using 
the well known background distribution, we calculate the probability that 
each event is a background signal, Bi = B{Ei), and the gaussian function 
Gi including the energy error SEi arising by the energy smearing. During the 
cooling phase, we assume that the muon and tau antineutrinos temperatures 
(T(P^) and T{i>r) respectively) are proportional with the electron antineutrino 
temperature (r(j?e)), i.e. T{vr)/T{ue) =T{v^,)/T{ve) = 1-2 [15]. 

At first neglecting the accretion phase, we solely consider the neutrino flux 
of cooling phase. We study the probability function L and we find a maximum 
when the model parameters reach the best-fit values shown in Table 1. We 



Table 1. Results for one-component (cooling) model with 2a errors 



Tc(MeV)|rc(sec)|_Rc(Km)|tK//(sec)|t/MB(sec)|tsafe(sec) 
^■^tVt |3.7+?1| 31+?^ I 0+°« I 0+«-^ I 0+4-^ 



remark that the best-fit value for Re is larger than the theoretically expected 
one (namely Re ~ Rns — lOKm [16]). We calculate the total energy carried 
by neutrinos during this phase corresponding to the gravitational binding 
energy of neutron star £{,. Using the equipartition hypothesis the relations 
£b = 6- £c{ve) = 3.39 IQ-^RlT^Tc = 3.87 • lO^^g^^ j^oy^ .^^j^gre the mean 
values of antineutrino energy are (Ep^) = IQMeV and (i?p^) ~ l2MeV, a bit 
lower than expected [15, 17]. 

Motivated by the experimental fact that about 40% of the SN1987A 
events have been recorded in the first second, we consider the accretion 
phase completing the emission model. We set Ma = 0.5Mq that is a rea- 
sonable value of the outer core mass, therefore we maximize the likelihood 
as a function of the other parameters. We find the best-fit values shown 



Table 2. Results for two-components (accretion and cooling) model with 2a errors 

Tc(MeV)|rc(sec)|i?c(Km)|ra(MeV)|ra(sec)|tif/i(sec)|t/Ms(sec)|tsafe(sec) 
r -1+2.1 \aa+3.6\ 10+18 I 91+0.2 In 7+1-3 1 n+0-** I n+"-7 I n+o-B 



in Table 2. Note that Tc, Re and Tq are very close to theoretical values 
expected [16]. The binding energy is the sum of two terms, the energy of 
neutrino emitted in the cooling phase, £c = 1.76 . lO^^erg, and the energy 
£a = 2 ■ £a{ve) = 4,.UMaT^Taip = 6.3 . lO^^^g^^ Carried by i^e and De in 



J ,j ^ , -K Neutrino Luminosity 

Log(L) (erg sec ) ' 

10"(51) 
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Neutrino mean energy 
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Fig. 1. The luminosity (a) and the mean energy (b) of Ve and Vx (dashed line and 
solid line, respectively) obtained from the emission model (accretion and cooling) 
exploiting the best fit values of Table 2. 
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Fig. 2. Time integrated energy spectra of neutrino for accretion (dashed line) and 
cooling (solid line) emission phase. 



the accretion phase, where (p = L dx cxp(— x^°)/(l + x Ta/0.5). Hence 
we obtain £b == 2.4 • lO^^g^^^ {Epja = lO.SMeV, {E^Jc = U.GMeV and 
(Ep^) — 15MeV. We report the neutrino luminosities, the neutrino energies 
mean values (as a function of time) in Fig. 1 (a) and (b), respectively. More- 
over, we plot the neutrino energy spectra in Fig. 2. 



4 The evidence for the phase of accretion 

In the analysis with accretion and cooling phases, we find that the absolute 
value of the likelihood function (in the best-fit point) is about 1000 times 
larger than the corresponding value in case of no accretion phase, giving a 
significant hint for an accretion phase. In fact, let us assume as null hypothesis 
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Hq the case where the accretion is absent and compare it with the alternative 
hypothesis Hi with an accretion phase described by two additional parameters 
Ta and r^. When we add 1/ = 2 degrees of freedom, we expect that the x^ 
will decrease by a certain amount Ax^- In order to determine the rejection 
interval for the hypothesis Hq, we perform a likelihood ratio test; then the 
required probability distribution function is the regularized gamma function 
Q{v/2,Ax^/2). When we go from Ma = (no accretion) to Ma = 0.5 M© 
(our reference point) the x^ diminishes by Ax^ — 13.4. As consequence, we 
reject the null hypothesis in favor of the hypothesis that accretion occurred 
with a significance of a = exp(— Z\x^/2) = 1.2 x 10""^. In Gaussian language, 
this amounts to 3.2(t. 
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